Quantum electrodynamic calculation of the hyperfine structure of ^He 
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I. INTRODUCTION 

Spectroscopic measurements of the helium atom have 
presently reached the level of accuracy at which they are sensi- 
tive to uncertainties of fundamental constants and the nuclear 
effects. This fact can be utilized for improving our knowledge 
on fundamental constants and nuclear parameters. A recent 
example is the independent determination of the fine structure 
constant a [1] made by comparison of the theoretical predic- 
tion and the experimental results for the ''He fine structure. 
Good agreement of the obtained value of a with the more pre- 
cise determinations [2, 3] provided a highly sensitive test of 
consistency of different theories across a wide range of en- 
ergy scales. 

Recent optical measurements in "^He [4, 5] and ^He [6] 
achieved the relative precision of about 10^^^ and thus be- 
came sensitive to the uncertainties of the Rydberg constant 
and the nuclear charge radius. These experiments created a 
possibility for the spectroscopic determination of the nuclear 
charge radii of -^He and '*He, with a significantly higher accu- 
racy than can be reached by the electron scattering methods 
[7]. Such determination is of particular interest now, in view 
of the discrepancy for the proton charge radius observed in 
the muonic hydrogen experiment [8] and the proposed follow- 
up experiment on the muonic helium [9]. Realization of this 
project requires progress in theory, namely a complete cal- 
culation of the m (J corrections to the energy levels. Such 
calculation is difficult but probably feasible, at least for the 
low-lying triplet states. 

While the present theory is not accurate enough to provide 
the nuclear charge radii of "^He and ^He separately, it can pro- 
vide their difference br^, as the isotope shift is considerably 
simpler to calculate than the energy levels. Two such deter- 
minations have recently been reported by experiments of the 
Amsterdam [4] and the Florence [5] groups, their results be- 
ing in disagreement of four standard deviations. There was 
also the older spectroscopic value of Sr^ by Shiner et al. [10], 
which relied on the theory of the hyperfine splitting available 
at that time. 

The main goal of the present investigation is to calculate the 
combined fine and hyperfine structure of the 2^P levels of ■^He 
with the complete treatment of the mo^ corrections. The re- 



sults obtained represent an order-of-magnitude improvement 
over the previous theory [11, 12] and are in good agreement 
with the experimental data available [5, 13]. The improved 
theory allows us to make a reevaluation of the Jr^ determi- 
nation by Shiner et al. [10], in reasonable agreement with the 
Florence result [5]. We also update the previous results for 
the hyperfine splitting of the '2?'S state of ■^He [14], making its 
treatment fully consistent with that of the !?¥ states. 



II. GENERAL APPROACH 

The !?¥ energy level in '^He is split by the hyperfine and 
fine structure effects. The hyperfine splitting is induced by 
the interaction between the dipole magnetic moment of the 
nucleus and that of the electrons, whereas the fine structure 
is due to the interaction between the electron spin and the 
electron orbital angular momentum. For heavy atoms, it is 
common that the hyperfine splitting (being suppressed by the 
electron-to-proton mass ratio) is much smaller than the fine 
structure splitting, and so the hyperfine and fine structure ef- 
fects can be investigated separately. For the 2'^P states of "^He, 
however, both effects are of the same order of magnitude and 
thus should be calculated together 

The general method of calculation of the combined fine and 
hyperfine structure is the nonrelativistic perturbation theory 
for quasidegenerate states. The active space of (strongly in- 
teracting) quasidegenerate states is defined as i^Pj , where F 
is the total angular momentum of the atomic state and J ~ Q, 
1, 2 is the electronic angular momentum. In the case of "^He, 
the nuclear spin / = ^/2, so i^ = J ± 1/2 > and the ac- 
tive space consists of five levels. Note that, unlike in previous 
studies [11, 12, 15], we do not include the 2^P state in the ac- 
tive space, as its mixing with the 2^P levels is relatively weak. 
The contribution of this state is accounted for by the standard 
perturbation theory for non-degenerate states. 

The energies of the 2^Pf states are the eigenvalues of the 
5x5 matrix of the effective Hamiltonian H, whose elements 
are 



EIj, = {FJMf\H\FJ'Mf) 



(1) 



where Mp is the projection of the total angular momentum 
F. (Since the energies do not depend on Mp, it can be fixed 



arbitrary.) In practical calculations, it is convenient to con- 
sider the shifts of the individual i^Pf levels with respect to 
the 2^P centroid. In other words, we require that 



Y,{2F+1)Y,eIj^Q. 



(2) 



In this case, all effects that do not depend on the nuclear and/or 
electron spin do not contribute and can be omitted in the cal- 
culations. 

The matrix elements of the effective Hamiltonian (1) are 
obtained in this work by expansion in terms of the fine struc- 
ture constant a. 
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where Hfs is the effective operator responsible for the fine- 
structure splitting in the absence of the nuclear spin and 
the other terms are the nuclear-spin-dependent contributions. 
-ffj^fg is the leading hyperfine Hamiltonian with the recoil 
and anomalous magnetic moment additions; it is of nominal 
order raa^ and contains parts of higher-order contributions. 
H^^l is the effective Hamiltonian of order ma^ and is derived 
in the present work. The next two terms in Eq. (3) are the 
second-order corrections that also contribute to order ma^. 
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ment additions, H^^ is the effective Hamiltonian responsible 

for the leading fine-structure splitting of order ma^, and H^^fl 
is the effective Hamiltonian responsible for spin independent 
effects of order ma^. Finally, -ffnuci&ho represents the nuclear 
effects and the higher-order [^ ma^) QED corrections pro- 
portional to the 5 function at the origin. This part cannot be 
accurately calculated at present, because of insufficient theo- 
retical knowledge of the nuclear structure. It will be obtained 
from the experiment on the ground-state hyperfine splitting in 
^He. 

In order to facilitate the evaluation of the matrix elements 
in Eq. (3), it is convenient to factorize out their dependence 
on the nuclear degrees of freedom F and Adp and on the elec- 
tronic angular momentum J. It can be achieved by observing 
that any operator contributing to Eq. (3) can be represented 
in terms of six basic angular-momentum operators: {S ■ L), 
(/ • L), (/ • 5), (5*5^)(2)(L*iJ')(2), PS^{UU)^^\ and 
P U {S"^ S^Y^\ where the second-order tensors are defined 
by {U LJ)(2) EE 1/2 L'L^ + 1/2 L^L' - 1/3 L^S'^ and the sum- 
mation over the repeated indices is assumed. Since the nuclear 
spin of helion / — 1/2, there are no operators quadratic in /. 
So, Eq. (3) can be represented as 



E 



J. J' 



+AM- S) + A{I- L) + A,u {P S^ {V L^)(2)) 
+A,,i(ri-'(5^5-')(2)), (4) 



where the constants Ai do not depend on F, Mp, J, J', and 
are represented by expectation values of purely electronic op- 
erators between the spatial 2^P wave functions. The matrix 
elements of the basic angular-momentum operators are calcu- 
lated analytically by means of the Racah algebra and listed in 
Appendix A. 

The first two terms in the right-hand-side of Eq. (4) do not 
depend on the nuclear spin and correspond to the fine struc- 
ture. The fine structure splitting in the absence of the nuclear 
spin was investigated in our previous investigations [1, 17], so 
we use results obtained there. 



III. LEADING HYPERFINE STRUCTURE 

The leading hyperfine-structure Hamiltonian H^^^ ' is well 
known. For our purposes, it is convenient to write it in the 
following form. 



7j(l+) = ^(l + «)a4 



'hfs 



mM 
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I-SaQa + P S\ Q'i 



(5) 



where S and Sa are the electron spin operators, S ~ {a\ + 
CT2)/2 and Sa — (f?! — '?2)/2, rrir is the reduced mass of 
the electron-nucleus system, M is the nuclear mass, and k is 
the magnetic moment anomaly related to the nuclear dipole 
magnetic moment /i by 
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where /is is the Bohr magneton. The electronic operators are 
given by (in atomic units) 



Q =(\^a,)-^^{b\T^)^^\v2)\, (7) 

Qa = (1 + fle) I 47r[(53(ri) - ^^(ra)] , (8) 



Q 



m 1 + 2k Z 



M l + K 2 



X P2 




r^i , ^2^ 


3 


„3 
^2. 



X(P1+P2), (9) 



Q' 



-(l+fle) 



z 



z 



S^3 



'I'l 



S'^ - 3 



2'2 



QJ =-(l + ae)^ 



(10) 



(11) 



where a^ is the anomalous magnetic moment of the electron. 



TABLE I: Expectation values of the ma hyperfine operators for the 
2 P state of He, with the mass polarization included, in a.u. The 
expected numerical uncertainty is less than 1 on the last significant 
digit. 
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21.092 193 
0.277 443 

-0.060 491 

0.140 325 



The matrix element of the Hamiltonian between the 2^Pj 



states can be represented as 



Derivation of the effective Hamiltonian to order ma" for an 
arbitrary state of helium is given in Appendix B. The result is 

^3 
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where 



P,^a^z{\n2-'^^A^[5\ri) + 5\r2)\ , (18) 
is the effective operator induced by the radiative QED effects, 
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where the shorthand notations for the matrix elements of the 
electronic operators are 

{^m) = (0.|Q|V»), (13) 

{$\Qm = -*6*-''"(0^iQj#fc), (14) 

(0IQI^) = {(f>Mj\i^j), (15) 

and the spatial triplet odd P wave function is represented as 



ip\ri,r2,r) = r\f{ri,r2,r) - (ri O r2) 



(16) 



with f{ri,r2,r) being a real scalar function of ri, r2, and 
r = \ri — f2\- The P-state wave function defined above is 
real and normalized to {P\P) — {Pi\Pi) — 1. 

Matrix elements of the operators Q, Q, and Q in Eq. (12) 
are evaluated between the wave functions that include the 
mass polarization term in the nonrelativistic Hamiltonian. Nu- 
merical results for the expectation values of these operators 
are presented in Table I. We observe that the Fermi contact 
interaction yields a dominating contribution, which could be 
aniticipated, as it comes from both Is and 2p electrons. This 
also explains why the hyperfine splitting is of the same order 
as the fine splitting, since the latter comes mainly from the 2p 
electron. 



IV. ma" CORRECTIONS 

Calculation of the ma^ contribution to the hyperfine struc- 
ture of the 2^P state is the principal objective of this work. 
Analogous calculation for the 2''5i state have already been 
performed in Ref. [14]. Here we verify the calculation for the 
2^5* state and extend it to the 2^P state. 
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It can be immediately seen that the operator Pnrad involves 
highly singular operators, whose matrix elements between the 
■^P functions diverge. However, the divergence cancels out if 
^nrad IS Considered together with the second-order ma^ con- 
tribution, as is explained below. 

The second-order ma^ correction SEscc is given by 



6Escc - 2 (H^f^ _ [H^f^ 



H^)^ 



(22) 



where -ff^fi 1^ obtained from H^^^ ' from Eq. (5) by dropping 
the electron magnetic moment anomaly a,, and the recoil part, 
and 



^(4) 
'nfs 



4'*^ = T + 5 • f + S' S^ T'^ +SA-fA, (23) 
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where H^^^ = T is the spin independent part of the effective 
Hamiltonian of order m,a*. 
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and the remaining operators are responsible for the fine struc- 
ture to order m,a^. 
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TABLE II: Expectation values of the ma*^ hyperfine operators for 
the 2 P state of He, in a.u. The expected numerical uncertainty is 
less than 1 on the last significant digit. 
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-^rad 

P' 

^ nrad 

P 
P 



-76.220 978 

0.367 674 

-0.529 385 

-0.361 070 



TABLE III: Individual second-order matrix elements, for different 
types of intermediate states "L, in a.u. The expected numerical un- 
certainty is less than 1 on the last significant digit, when not given 
explicitly. 
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Let us now consider the sum of singular contributions. 
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While the expectation values of the operators Q and T are 
finite, the second-order matrix element of these operators di- 
verges. In order to eliminate the divergences, we regularize 
the Coulomb electron-nucleus interaction by introducing an 
artificial parameter A, 



Z Z , 
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(29) 



All other electron-nucleus interaction terms are regularized in 
the same way. This entails the following replacements in ef- 
fective Hamiltonians, 



AnZS'^ira) = - V' — ^ -V' — (l - e-^"^"" 
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Once the electron-nucleus interaction is regularized, one can 
in principle calculate all matrix elements and take the limit 
A — > oo. However, since matrix elements can be calculated 
only numerically, it is more convenient to transform the ef- 
fective operators to the regular form, where A can be taken to 
infinity before numerical calculations. 

To this end, we transform the operators in the second-order 
matrix element T ^ T' and Q ^ Q' by 



t'^-t-\i:{^^^e-h]. 



Q' ^ Q + \Y.[^^E-H 



(32) 
(33) 



where {...,...} is the commutator and the implicit A- 
regularization is assumed. After this transformation, the sin- 
gular part takes the form 



?i ■ E ■ 



mM 



(l + K)a'^I-S 



I pi \ 

\ nrad/ 



1 



{E-H)' 



T' + h.c. 



(34) 



where both the first and second order terms are separately fi- 
nite and thus the limit A ^ oo can be evaluated analytically. 
The result for the regularized operator P^-^^^^ is 
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Matrix elements of the regularized ma^ Hamiltonian are ob- 
tained by 
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The numerical results for the expectation values of the elec- 
tronic operators are listed in Table II. 

The second-order ma^ correction [Eq. (22)] is conceptually 
simple but its calculation is rather involved technically. This is 
partly because of the coupling with intermediate states of dif- 
ferent symmetries {^P, ^P, ^D, ^D, and "^F) and partly due 
to the presence of nearly singular operators. For achieving nu- 
merically stable results, some operators had to be transformed 
to the regular form by using the same method as for the sin- 
gular part. Details of the calculation and the regularization 
procedure are described in Appendix C. 

Numerical results for the second-order ma^ corrections 



are presented in Table III. Note a large contribution of the 
^P intermediate states, which is mainly due to the 2^P-2^P 
mixing. We also mention that the numerical uncertainties of 
the second-order corrections are completely negligible on the 
level of the total theoretical error. 



V. SECOND-ORDER HFS CORRECTION 

The fifth term in Eq. (3) is the second-order hyperfine cor- 
rection. Its nominal order is {m? /M) a^ . However, it is en- 
hanced by the small 2^P-2^P energy difference in the de- 
nominator, which makes it numerically significant. Rigorous 
calculation of this contribution is difficult because of diver- 
gences, which are canceled by the corresponding contribution 
from the nuclear-structure effects. In the present work, we ap- 
proximate the second-order hyperfine contribuiton by keeping 
only the lowest lying ^P/^ intermediate state in the sum over 
spectrum. 

The second-order hyperfine correction induces contribu- 
tions to the isotope shift (i.e., to the centroid of the 2^P level). 
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and to the hyperfine structure. 
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TABLE IV: Second-order hyperfine matrix elements, in a.u., 
all figures shown are exact. Notations are: (A, B) = 

{2-^ P\A\2^ P) {2^P\B\2-^P) /[E{2^P) - E{2^P)]. The expected 
numerical uncertainty is less than 1 on the last significant digit. 



(Qa , Qa) 


-47 774.980 


(Qa , Qa) 


249.528 


{Qa , Qa) 


-1.303 



Numerical results for the second-order hyperfine corrections 
are presented in Table IV. 

The approximate treatment of the second-order hyperfine 
contribution yields, in our opinion, the dominant uncertainty 
of our theoretical predictions for the hyperfine structure. It 
is, therefore, important to estimate the neglected part. To this 
end, we introduce a different approximation for the second- 
order hyperfine contribution, which is obtained from Eq. (39) 
by the following substitution 
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where Q^ is the regularized 5 function operator given by 
Eq. (C20). The difference between these two approximations 
(of about a half percent) is used as the estimated error of the 
theoretical hyperfine structure. 



VI. NUCLEAR-STRUCTURE CONTRIBUTION 



In order to infer the nuclear-structure contribution, we 
subtract from the ■^He+ experimental result the hydrogenic 
limit of all corrections accounted for in the previous sections. 
The remainder is the nuclear-structure contribution plus some 
(small) higher-order QED corrections. Specifically, we de- 
fine the nuclear-structure and higher-order QED contribution 

^nucl&ho ^S 
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m M 
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2 \ 2 
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(41) 



where E^^p{ls,^ Hc+)//i == 8665649.867 (10) kHz is the ex- 
perimental result of Ref. [16]. The above equation yields 



a 



nucl&ho 



-0.031891 



(42) 



This value does not have any theoretical uncertainty (by def- 
inition) and thus is accurate to all figures given. The corre- 
sponding contribution to the hyperfine structure in neutral he- 
lium then is 
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The last term in Eq. (3) is the nuclear-structure and higher- 
order QED contribution. Accurate calculation of the nu- 
clear contribution is presently not possible due to insufficient 
knowledge of the nuclear dynamics. However, one may claim 
that the dominant part of this effect comes from a local op- 
erator proportional to the S function at the origin. Using 
this assumption, one can infer the nuclear-structure correc- 
tion in neutral ■^He from the experimental value of the ground- 
state hyperfine splitting in ■^He+, measured very accurately by 
Schliisser etal [16] half a century ago. This approach has been 
used also in the previous studies of the "^He hyperfine structure 
[11,12]. 



VII. SMALL CORRECTIONS 

Here we pick up some higher-order contributions, which 
are not accounted for by Cnuci&ho but might be relevant for 
comparison with experiment. Since the mixing of the 2^P and 
2^P levels by the fine-structure operator is large, we calculate 
here the anomalous magnetic moment and recoil corrections 
to this mixing, which are of orders a^ w? /M and a® w? /M"^, 
respectively. 

The nuclear spin dependent 2'^P-2^P fine- structure mixing 
is given by 
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(44) 



where the Qa and Ta operators include the anomalous mag- netic moment and the recoil additions, 

QA = {l + ae)^47r[d^{n)-S^{r2)] 



Ta = {1 + 2a, 
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TABLE V: Nuclear-spin dependent fine-structure mixing matrix ele- 
ments, with the anomalous magnetic moment and mass polarization 
(upper entry) and without (lower entry), in a.u., all figures shown 
are exact, (A, B) = {2-^P\A\2^P) {2^P\B\2^P) /[E{2^P) - 
E{2 P)]. The expected numerical uncertainty is less than 1 on the 
last significant digit. 



{Qa , Ta) 
{Qa , Ta) 



156.1127 
155.9034 
-0.81535 
-0.81428 



and the matrix elements are calulated between the wave func- 
tions that include the mass polarization correction. 

The higher-order mixing contribution is obtained from 
Eq. (44) after subtracting the part that is already accounted for 
[namely, the n ~ 2 term in Eq. (C6)]. Numerical results for 
the corresponding matrix elements are given in Table V. The 
recoil correction is due to both the the mass polarization and 
the recoil addition to Ta- Further corrections to this mixing 
(e.g., those coming from higher-order relativistic corrections) 
are not known and contribute to the uncertainty of final results. 



TABLE VL Theoretical results for individual 2^Pf levels in *He, 
relative to the 2^P centroid energy, in comparison with available ex- 
perimental data, in kHz. The first error of the theoretical values is the 
uncertainty due to the hyperfine structure and the second error is the 
uncertainty due to the fine structure. 



3r,F=l/2 



2-'P, 



3 pF=3/2 



2-'Pr 



23pF=l/2 



3pF=3/2 



2^Pl 



3pF=5/2 



2^P, 



27923 393.7 
27923 394.7 
27923 398.3 

498 547.3 
498 543.7 
498 547.3 

-169 462.2 
-169 463.3 
-169 460.2 

-4 681676.3 
-4 681676.3 
-4 681672.1 

-6 462 557.8 
-6 462 555.3 
-6 462 562.8 



(0.2) 


(2.2) 


(1.9) 


(1.4) 


(2.1) 


(2.1) 


(0.5) 


(1.7) 


(1.8) 


(0.7) 


(1.5) 


(1.6) 


(0.7) 


(1.5) 


(1.6) 



(2.5) 



(0.4) 



(0.8) 



(0.2) 



(1.0) 



Smiciklas [13] 
Cancio et al. [5] 



Smiciklas [13] 
Cancio et al. [5] 



Smiciklas [13] 
Cancio et al. [5] 



Smiciklas [13] 
Cancio et al. [5] 



Smiciklas [13] 
Cancio et al. [5] 



VIIL RESULTS AND DISCUSSION 



where 



In this section we present the total theoretical predictions 
for the mixed fine and hyperfine structure of the 2^P states in 
•^He. The numerical values of nuclear parameters used in the 
calculations are [18] 

fi/fiB = -1.158 740 958(14) 10"^ (47) 

m/M = 1.819 543 0761(17)10"'^. (48) 

The conversion factors relevant for this work are 



to: 



toM 



mM 
m 



f7(l + «)"' 



-202 887.3247 kHz x h, (49) 



{1 + K)a^ = -10.8040 kHz X /i, 



2 „,6 



m? M"^ 



17^(1 + -r" 



0.0063 kHz X /i , 



(50) 



(51) 



where h is the Planck constant. 

The fine structure splitting in the absence of the nuclear 
spin was calculated in our previous investigations [1, 17], with 
numerical results reported for ^He. In this work, we reevaluate 
all nuclear-mass-dependent corrections to the fine structure, in 
order to extend our calculation to ■^He. The numerical results 



(^f.),/=o =9(8./oi + 5/i2), 

(^ft),7=l =^(-./01+5/l2): 



(52) 
(53) 
(54) 



/oi ==29 616 676.5 (1.7) kHz, 
/i2== 2 292 167.6 (1.7) kHz. 



(55) 
(56) 



We now have all contributions to the elements of the Hamil- 
tonian matrix Ej j, in Eq. (3). Diagonalizing the matrix, we 
obtain the positions of the energy levels of the 2^Pf states in 
■^He, relative to the centroid of the 2^P level. The numerical 
results are listed in Table VI for the individual energy levels 
and in Table VII for the transitions between the fine and hyper- 
fine levels. Our theoretical values have two uncertainties, the 
first one being due to the hyperfine effects and the second one, 
due to the fine-structure effects. The nuclear-spin dependent 
effects are calculated with an accuracy of about 1 kHz (and 
even better in some cases). This accuracy is limited mainly by 
the incomplete treatment of the second-order hyperfine cor- 
rection. The second uncertainty of the theoretical energies is 
due to the fine-structure effects. It comes from the 1.7 kHz 
error of /oi and /12 in Eqs. (55) and (56), which is exactly 
the same as for the fine structure in ''He. Interestingly, the 
sensitivity of different levels to the error of the fine-structure 
effects is very much different, varying from 0.2 kHz for the 



,3/2 



1/2 



2*P'{' level to 2.5 kHz for 2'^P^' . As could be expected, 
the transitions between the levels with the same value of J are 
less sensitive to the error of the fine-structure effects than the 
J — J' transitions. 

The present theoretical results can be compared with the 
experiment by Smiciklas [13] and with our recently reported 
absolute frequency measurements of the 2'^S'-2'^P transitions 
[5]. The latter experiment was carried out by using the op- 
tical frequency comb assisted multi-resonant precision spec- 
troscopy [19] and by measuring simultaneously both optical 



TABLE VII: Experimental and theoretical hyperfine transitions, in 
kHz. The first error in the theoretical prediction is the uncertainty 
due to the hyperfine splitting and the second error is the uncertainty 
due to the fine-structure splitting. 



{J,F)-{J',F' 



Value 



Reference 



(0,1/2) -(1,1/2) 

(0,1/2) -(2, 3/2) 
(2, 3/2) -(1,1/2) 

(1,3/2) -(2,5/2) 

(1,1/2) -(1,3/2) 



28 092 855.9 (0.5) (1.7) 

28 092 870 (60) Morton et al. [11] 

28 092 858 (3) Smiciklas [13] 

28 092 858.6(2.1)" Cancio et al. [5] 

27424 846.4 (1.4) (2.9) 

27424 837 (12) Wu&Drake [12] 

27424 851.0(3.0) Cancio e? a/. [5] 

668 009.5 (1.4) (1.2) 

668 033 (9) Wu&Drake [12] 

668 007(3) Smiciklas [13] 

668 007.5(3.2) Cancio ef a/. [5] 

1780 881.5 (1.0) (1.2) 

1 780 880 ( 1 ) Wu&Drake [12] 

1780 879(3) Smiciklas [13] 

1780 890.7(3.5) Cancio ef a/. [5] 

4512 214.1 (0.8) (0.5) 

4 512 191 (12) Wu&Drake [12] 

4 512 213(3) Smiciklas [13] 

4 512 211.9 (2.7) Cancio et al. [5] 



(2, 3/2) - (2, 5/2) 6 961 105.1 (1.5) (0.5) 

"Averaged value of the two differences between the measured optical tran- 
sitions. 



and microwave hyperfine transition frequencies. Agreement 
found between the microwave frequencies and the difference 
of the optical transition frequencies was used as a confirma- 
tion of the obtained experimental results. Comparison of the 
two independent measurements (see Table VII) shows good 
agreement in all cases except for the P{ -Pj transition. 

The experimental values for individual energy levels listed 
in Table VI are obtained from the transition energies reported 
in the original references [5, 13] by using the definition of the 
centroid [Eq. (2)] and the experimental hyperfine shift of the 
2'^S state [16]. We observe that theoretical and experimen- 
tal results are at the same level of accuracy of about 2-3 kHz 
and in very good agreement with each other. It can be con- 
cluded that our calculation of the ma^ correction represented 
an important advance over the previous theory [11, 12] and 
significantly improved agreement between theory and the ex- 
periment. 

Table VIII shows our numerical results for the hyperfine 
splitting of the 2'^S state of '^He. The results listed represent 
an update of the calculation described in detail in Ref. [14]. 
As compared to that work, we have (i) slightly improved the 
numerical accuracy and (ii) made the treatment of the second- 
order hyperfine and the nuclear-structure contributions to be 
fully consistent with that for the 2^P states. The uncertainty 
of the theoretical prediction comes from the second-order hy- 



2'S-2'S IS [4,5] . 



2=P-2'S IS [10 and this work] 
"" 2'P-2'S IS [5] 

Electron-scattering [7,25] 



Nuclear theory [23,24] ► 



0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 

5r'['He-''He](fm') 

FIG. 1 : (Color online) Different determinations of the difference of 
the squared nuclear charge radii of '^He and ^He. 



TABLE VIII: Hyperfine splitting of the 2^S state of ^He, in kHz. 



/rr(4) „(4) „(4)^ 

l--"tifs'-"nfs + -"fs ) 

-^nucl&ho 

Total 
Experiment [20] 



6 740 451.46 
-1313.99 
2189.81 
-60.52 ±1.7 
-1566.83 
6 739 699.93 ±1.7 
6 739 701.177(16) 



perfine correction and was estimated in the same way as for 
the 2^P state. Very good agreement of the 2'^S theoretical re- 
sult with the experimental value [16] gives us additional con- 
fidence in our estimation of errors for the 2^P state. 

Having calculated the hyperfine and fine structure of the 
2'^P levels, we are now in a position to obtain an improved 
determination of "^He-^He nuclear charge radii difference Sr'^ 
from the isotope shift measurement by Shiner et al. [10]. 
In order to extract Sr^ from the measured energy differ- 
ence, £;(3Hc,23Pg/^-235'f/^)-£^(4Hc,23p2-23S'i), we sub- 
tract the experimental hyperfine shift of the 2'^S state [16], 

q 1/2 

the theoretical shift of the 2'^Pg ' level with respect to the 
centroid energy (obtained in this work), the theoretical fine 
shift of the 2^P2 level with respect to the centroid [1], and the 
theoretical isotope shift of the centroids for the point nucleus 
[5] (see Table IX). The remainder SE comes from the finite 
nuclear size effect and is proportional to the difference of the 
mean square charge radii, 6E — C Sr^, where coefficient C is 
evaluated in this work to be C = -1212.2(1) kHz/fm^. The 
resulting value 



5r^ = r2(^He) - r'^{'^YLe) = 1.066(4) W 



(57) 



is in reasonable agreement with the result by Cancio et al. [5], 
5r^ ~ 1.074(3) fm^, and in significant disagreement with the 
result of Ref. [4] (updated in Ref. [5] by recalculation of the 
isotope shift), 5r^ = 1.028(11) fm'^. Figure 1 shows graph- 
ically the comparison of different determinations of Sr^, in- 
cluding the results from nuclear theory [23, 24] and from nu- 
clear electron-scattering measurements [7, 25]. 



TABLE IX: Determination of the difference of the mean square nuclear charge radii Sr^ of '^He and ^He from the measurement by Shiner et 
al. [10]. The remainder SE is proportinal to 5r , 5E = C 5r , with C — —1212.2(1) kHz/fm , see text for details. Units are kHz. 



S(^He, 2^P^^^ - 2^Sf'^) - _B(^He, 2^P2 

5Shfs(23sf/') 

5EW2^P2) 

-5Shfe(23Po'") 

-SEiso{2^P - 2^S) (point nucleus) 

6E 



-2^Si) 



810 599.(3.) 

-2 246 567.059(5) 

-4 309 074.2(1.7) 

-27923 393.7(1.7) 

33 667143.2(3.9) 

-1292.8(5.2) 



Experiment [10] 

Experiment [16] 

Theory, [1] and this work, Eq. (54) 

Theory, this work. Table VI 

Theory [5] 



IX. SUMMARY 

In summary, we have calculated the mixed hyperfine and 
fine structure of the 2^P states of "^He. Our investigation ad- 
vances the previous theory by a complete calculation of the 
ma^ correction, which leads to an order-of-magnitude im- 
provement in accuracy. Theoretical predictions for most of 

Q 1/2 

the transitions are accurate to better than 2 kHz. The 2'^P^ ' - 
2^P^ transition is calculated up to 0.9 kHz, which is cur- 
rently the most precise theoretical result for the helium tran- 
sitions. Both the hypefine and fine-structure transitions are in 
good agreement with measurements by Cancio et al. [5] and 
by Smiciklas [13]. 

Since the present theoretical accuracy for the fine-structure 
transitions in ^He is comparable to that for ''He, one can 
in principle use the ■^He spectroscopy for the determination 
of the fine structure constant a, as in was done for "'He in 
Ref. [1]. However, such determination does not bring much 



advantage at present, since the experimental accuracy for ■^He 
is lower than for ''He [21, 22]. Because of this, we do not 
determine a from the ^He transitions in this work. However, 
if the experimental precision of the "^He hyperfine structure is 
improved to the sub-kHz level, one will be able to use these 
results for the spectropic determination of the fine structure 
constant. 



Another application of the hyperfine structure measure- 
ments in "^He might be determination of the dipole magnetic 
moment of helion. The principal problem here is that the 
present theory obtains the nuclear-structure contribution from 
the experimental value of the hyperfine splitting in '^He+. This 
greatly reduces the sensitivity of the final theoretical predic- 
tion on the nuclear magnetic moment. Our calculation shows 
that at present, the hyperfine structure of '^He allows for de- 
termination of the dipole magnetic moment of helion with an 
accuracy of about 5 x 10^^ only. 
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Appendix A: Angular-momentum matrix elements 



In this section we calculate the matrix elements of the basic angular-momentum operators that are relevant for the present 
work. The angular-momentum algebra is conveniently done with formulas from Ref. [26]. The results are 



{FIJLSMf\I-S\FIJ' 



'l's'Mf) = (-i)^+'^+^ { J' / ^ } ^^W+WT+T) 



{FIJLSMf\I ■L\FIJ 



'L'S'Mp) = {-iy+''+^ { J' / T } vW+WT+T) 



X Sl,l' Ss,s' (-1)^+5+'^'+i v/i(L + l)(2L + l)(2J+l)(2J' + l) \ j, f ^ 



(A2) 



{P S' (V i-'')(2)) = 1 [J{J + 1) + J'{.r + 1) - 2L{L + 1) - 2S{S + 1)] (f ■ i) - i L{L + 1) (/ • S) , (A3) 
{r L^ (5* 5^)(2)) = i [J{J + 1) + J'( J' + 1) - 2L(i + 1) - 2S{S + 1)] (f ■ S) -]- S{S + 1) (/ • L) , (A4) 






(A5) 



Appendix B: Derivation of H^ 



(6) 



We start with the Breit Hamiltonian Hbp of the atomic system in the external magnetic field, 



HbP — 2^ Ha + 2_^ ^ab 



a a.b:a>b 

2 TO Ta 2 ?Ti 8 m'^ 2 TO. 4 TO,"' r^ 

■2 I =;2 



8 m' 



(ffa • -Bq 7?^ + 7?^ (Tq • Sa 



-Hab = \ 7 [rab) - TT: 

rab m^ 
a 



' ab'ab 1 

^2 ) ^b 

' ab 



+ -. 5-^r i^a ■ rab X (2 TTb - TTq) - CTfc • fab X (2 7?^ - TTft) 

4 m'^r'lfj 

ao \ ab / 

where if = p — e A. Magnetic fields A and B induced by the nuclear magnetic moment are 

r 
47r '" r^ 



r,^ e^ r „ (1 + k) ^ r 

eA{r) = -— fix — = -Za ., ' I x 



M 



^S-irn = (vxl)-^-z.<l±^f .V,r + z„<l±^l(.«-3:lf^) 



p 



(Bl) 



(B2) 



(B3) 

(B4) 
(B5) 



The leading-order interaction between the nuclear spin / and the electron spin <?„ is obtained from the nonrelativistic terms 
1/2 m TT^ and e/2 m (Ta ■ Ba, yielding 



H, 



(4) 
hfs 



E 



6 ^ g -. 

— Pa • ^(^i) - 7; ^a ■ B{ra) 

TO 2m 



(B6) 
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in agreement with Eq. (5). The relativistic correction to the hyperfine interaction is similarly obtained from the relativistic terms 
in the Breit Hamiltonian Hgp, 



H, 



(6) _ 
hfs ~ 



Z a 



. r, 
a 



E i;;;! ^-^ • i >^ [-^^(^^)] + ^;3 (<?- • BaP^+f! Sa ■ b,) 



E 






8 m^ 
3a ■ rah X [-2eA{rb) + eA{ra)] 



4?Ti" 



Y: [Pa' Pa ■ A{ra) + Pa " Mr a) Pa'] " E ^ ^U £ + ^^ ) [- e A^ (n)] 



a^b 



' ab 



Using A and B from Eqs. (B4) and (B5), we derive Eq. (17). 



Appendix C: Second order matrix elements 



(B7) 
(B8) 



The second-order ma® correction 6Escc [see Eq. (22)] is split into several parts in accordance with the symmetry of the 
intermediate states. 



(^-Esnr — 



mM 



{1 + K)a^ 



SEsingCP) + SE^^g^P) + SECP) + SE^D) + 6E{^D) + 5E{^F) 



(CI) 



The angular momentum algebra for different intermediate states is simple but rather tedious and is performed with help of the 
symbolic algebra computer program. Below we list the resulting formulas expressed in a form convenient for the numerical 
evaluation. 

The singular part with the ^P intermediate states is given by (after removing all divergencies) 



<5£^si„g('P) = 



1 



(E-H) 



-T')2{I.S) 



where E = E{2^P), 



«'-<?+ E^.^-'^-^B-v. + S-v 



2Z fn 



(C2) 



(C3) 



1 V- f ^ 1 1 / Z Z 1 
T = T--y \—,E-H\ = -- (E+ — + 



P~^2-r^''' 



r* rJ 



Pi 






4 rS 2 r'^ 



(C4) 



and it is assumed that the operators T' and Q' act on the function on the right hand side that satisfies the Schrodinger equation 
with the energy E. 

The regular part with the '^P intermediate states is 



sE,,,ep) ^ j2 



n>2 



E{2^P) - E(v?P) 



{2^P\Q\n^P) {n^P\f\2^P) {- I ■ L + P U {S' S^)^^'>) 



+ {2^P\Q\n^P) {n^P\f\2^P) 
+ {2^P\Q\n^P) {n^P\f\2^P} 
+ {2^P\Q\n^P} {n^P\f\2^P) 
+ {2^P\Q\n^P} {n^P\T\2^P} 
+ (2^P\Q\n^P) {n^P\f\2^P) 
+ {2^P\Q\n^P) (n^P\f\2^P) 



■ P S^ {V L^)(2)) + {2^P\Q\n^P) {n^P\T\2'^P) {I- L) 



\i-s + \rsHL'uY^^) 



10 

6 



PU {S'S^f'^^) 
PS^VL^Y^^) 
I1.L + —P S^ {V i^')(2) -—P L^ {S' S-^y^)) 



1 ^ 91 97 

-I-S PS^ (V iJ)(2) -p V [S' S'J')(2)) 

5 100 ^ ^ 100 ^ ^ ' 



(C5) 
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The other parts are 



21) 2U 

<5ii;(3i?) = E E{2^P) - E{n^D) [^^'^'Q'"'^) (n'^^lf |2^P) {\ I-S-\ P ffl (^ L^r^) 

+ {2^P\Q\n^D) {n^D\f\2^P) {^ I ■ S + ^I' S^ {V L^)^^'> - -^ P L^ (S"' S')^^'>) 

9 30 10 

+ {2^P\Q\n^D) {n^D\f\2^P) (-^ / . £ - A /» s' {L' L^)(2) - ^ P L^' (5* S^)*^)) 

o J-U J-U 



(C6) 



(2^P\Q\n^D)(n^D\T\2^P) ( — P LUS' S^) 

5 



(2) 



(C7) 



5ECD) = E E{2^P)-E{n^D) (^'P\Q^\^'^) {n' D\fA\2' P) {-^ I ■ L-^P S^ {V L^^^^ + ^P U {S^ S^^^^) 



(C8) 



SECF) = E ^r93p^%r.3F^ (2^P|g|n3F) (n3p|f|23p)(i / • ^ - ^ P S^ {V L^f^ + ^ P U {S^ S^)^'^) 



E{2^P) - E{n^F) 



The matrix elements with the D and P-state wave functions are defined by 



10 



(</'IQI^) 



— jjk I 



{(t>^\QJl\1pkl) = e'^^MQiMi) = {i'\Q\4>) , 



((/)|(5|VJ) = {(t)i\Qjk\i^i;jk) = {i^ijk\Qij\4'k) = (^IQIf^) , 

where tj) denotes the odd D wave function, 
and V' denotes the odd F wave function, 



(C9) 

(CIO) 
(Cll) 
(C12) 

(CI 3) 



i, 



ijk 



r\r{ri - y {5^Jr\ + 5ikr\ + <5jferi) 



f{ri,r2,r) + - 



ijk, ijk, ijk 



%(r^,^ + 2r[r{r^) ^-^{r{r^, + 2r[r{ri) - ^-^(r\r\ + 2r{r{r\] 



g{ri,r2,r) ± (ri o ra) . 



(C14) 



All wave functions are real and normalized by (f/'^IV'*) = (V'*''IV'*"') = (^''■''^IV'*"''^) = 1- 

The second-order corrections listed above are finite. However, some matrix elements are too singular for the direct numerical 
evaluation and need to be transformed to a more regular form. The regularization method is described in Sec. IV. So, the 
operator P is transformed by Eq. (32), thus yielding 



E 



1 



n>2 



E{23P) - E{n^P 



■{2^P\Q\n^P){n^P\P\2^P) 



(C15) 

We observe that, after the transformation, the singular part of the operator P is absorbed in the first-order matrix element, 
whereas the second-order correction contains only the (more regular) operator T'. 
Similarly, the operator Q*-' is transformed by 



Q 



nj 



>- J Q .L^\ r3 r3 r3 r^ J 



(C16) 
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where E = E{2^P), 

1 Z / r^ r^ \ 

6Q'i = Y — \S'^ - 3^V^ (C17) 



6 ■'^-^ Ta V ^i 

and it is assumed that the function on the right hand side of Q' satisfies the Schrodinger ahgn with energy E. The second-order 
matrix elements with Q*-' then are transformed by 



y^ ^^/23p\Q\n3p\ {n^p\f\2^p) 



n>2 



E{23P) - E{n^P) 



= E J^(23P) - g(n3p) ^^'^'^''"'^^ (n3p|f |23p) + (-je'^'") (23p,|<5g,fc T, |23p„) - {2^P\5Q\2^P) {2^P\f\2^P) , 

(C18) 
and 

y^ "" — ^^l2ip\Q\n^p\ {n^p\T\2^P) = V , , , '^ — -^-- (2^ P\Q' \n^ P) {n^P\T'\2^P) 

^^E{2^P)-E{n^Py '^' '^ ' ' ' ^^E{2^P)-E{n^Py '^ ' '^ ^ ' ' ' 

\ z \ z 

+ {2^P\- ^ — Q\2^P) - {2^P\- ^ —\2^ P) {2^ P\Q\2^ P) + {2^P\SQT'\2^P) - {2^P\dQ\2^P) {2^P\T'\2^P) . 

(C19) 
The operator Qa is transformed as 



_ ^ , 2 I Z ^ „ ,,1 2Z / ri ^ , fa 



g:.^g. + 3|---,i^-i/|^-^-^-v, + -|.v.j, (C20) 

where it is assumed that the function on the right hand side of Q^ satisfies the Schrodinger equation with the energy E. The 
regularized form of the second-order matrix element then is 



y^ ^ ^ ^ ^ — -__ (23p|g l^ip) (niplf 123/3) 
^ E{2^P)-E{n^P)^ '^ ' '^ ' ' ' 



E Ei2^P)-Ein^P) ^^'P\Q'-\-'P) (-'P\TA^'P) - I (2^^1 (f - f ) TA^''^) • (C21) 

I 



